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The  receptivity  problem  for  Gortler  vortices  induced  by  wall  roughness  is  investigated. 
The  roughness  is  modelled  by  small  amplitude  perturbations  to  the  curved  wall  over  which 
the  flow  takes  place.  The  amplitude  of  these  perturbations  is  taken  to  be  sufficiently  small  for 
the  induced  Gortler  vortices  to  be  described  by  linear  theory.  The  roughness  is  assumed  to 
vary  in  the  spanwise  direction  on  the  boundary  layer  lengthscale,  whilst  in  the  flow  direction 
the  corresponding  variation  is  on  the  lengthscale  over  which  the  wall  curvature  varies.  In  fact 
the  latter  condition  can  be  relaxed  to  allow  for  a  faster  streamwise  roughness  variation  so  long 
as  the  variation  does  not  become  as  fast  as  that  in  the  spanwise  direction.  The  function  which 
describes  the  roughness  is  assumed  to  be  such  that  its  spanwise  and  streamwise  dependences 
can  be  separated;  this  enables  us  to  make  progress  by  taking  Fourier  or  Laplace  transforms 
where  appropriate.  The  cases  of  isolated  and  distributed  roughness  elements  are  investigated 
eind  the  coupling  coefficient  which  relates  the  amplitude  of  the  forcing  and  the  induced 
vortex  amplitude  is  found  asymptotically  in  the  small  wavelength  limit.  It  is  shown  that 
this  coefficient  is  exponentially  small  in  the  latter  limit  so  that  it  is  unlikely  that  this  mode 
can  be  stimulated  directly  by  wall  roughness.  The  situation  at  0(1)  wavelengths  is  quite 
different  and  this  is  investigated  numerically  for  different  forcing  functions.  It  is  found  that 
an  isolated  roughness  element  induces  a  vortex  field  which  grows  within  a  wedge  at  a  finite 
distance  downstream  of  the  element.  However,  immediately  downstream  of  the  obstacle  the 
disturbed  flow  produced  by  the  element  decays  in  amplitude.  The  receptivity  problem  at 
larger  Gortler  numbers  appropriate  to  relatively  large  wall  curvature  is  discussed  in  detail. 
It  IS  found  that  the  fastest  growing  linear  mode  of  the  Gortler  instability  equations  has 
wavenumber  proportional  to  the  one  fifth  power  of  the  Gortler  number.  The  mode  can  be 
related  to  both  inviscid  disturbances  and  the  disturbances  appropriate  to  the  right  hand 
branch  of  neutral  curve  for  Gortler  vortices.  The  coupling  co-efficient  between  this,  the 
fastest  growing  vortex,  and  the  forcing  function  is  found  in  closed  form. 
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No.  NASl-18605  while  the  authors  were  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. 
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Introduction 

Our  concern  is  with  the  mechanism  by  which  Gortler  vortices  in  incompressible  bound¬ 
ary  layers  are  stimulated  by  wall  roughness  elements.  We  shall  consider  the  cases  of  isolated 
and  distributed  roughness  over  the  whole  range  of  unstable  Gortler  numbers.  In  a  recent 
paper,  Hall  (1990),  the  Gortler  receptivity  problem  for  free-stream  disturbances  was  dis¬ 
cussed;  the  present  paper  completes  the  discussion  of  the  linear  receptivity  problem  for 
Gortler  vortices.  However,  some  of  the  cases  considered  in  this  paper  can  be  generalized 
to  take  account  of  finite  amplitude  effects  but  those  situations  will  not  be  discussed  here. 
Before  discussing  the  receptivity  problem  in  more  detail  we  shall  remind  the  reader  of  the 
main  results  of  linear  stability  theory  applied  to  boundary  layers  on  curved  walls.  The 
paper  mentioned  above.  Hall  (1990),  is  for  the  most  part  a  review  of  the  linear  and  non¬ 
linear  stages  of  Gortler  vortices  so  the  reader  is  referred  to  that  paper  for  a  more  detailed 
discussion  of  the  various  stages  of  Gortler  vortex  growth. 

The  main  feature  of  Gortler  vortices  which  makes  them  behave  in  a  quite  different 
manner  than  for  example  Tollmien-Schlichting  waves  or  Rayleigh  waves  is  that  they  are 
almost  always  dominated  by  nonparallel  effects.  This  means  that,  at  0(1)  wavenumbers 
and  Gortler  numbers,  the  concept  of  a  unique  neutral  curve  is  not  tenable;  the  position 
where  a  given  vortex  begins  to  grow  is  a  function  of  its  upstream  history.  Mathematically 
this  property  manifests  itself  through  the  parabolic  nature  of  the  disturbance  equations  for 
Gortler  vortices.  Thus  Hall  (1983)  showed  that  at  order  one  values  of  the  Gortler  number 
and  wavenumber  a  small  amplitude  Gortler  vortex  is  described  by  a  parabolic  system  of 
equations  which  in  general  must  be  solved  numerically  for  each  vortex  wavenumber.  The 
neutral  curve  corresponding  to  any  fixed  initial  disturbance  can  be  computed  by  marching 
the  disturbance  equations  downstream,  but  it  is  a  function  of  the  initial  disturbance. 

The  only  exception  to  the  situation  discussed  above  is  that  when  the  vortex  has  a 
wavelength  small  compared  to  the  boundary  layer  thickness.  In  this  regime  the  vortex 
‘feels’  the  local  structure  of  the  boundary  layer  and  is  therefore  able  to  develop  in  a  quasi¬ 
parallel  manner.  Hall  (1982a)  showed  that  it  is  then  possible  to  define  a  unique  right  hand 
branch  of  the  neutral  curve  which,  at  zeroth  order,  has  the  Gortler  number  proportional 
to  the  fourth  power  of  the  vortex  wavenumber. 

It  had  been  assumed  by  authors  previous  to  Hall  (1982a)  that  Gortler  vortices  could 
be  described  in  a  self-consistent  manner  by  making  the  parallel  flow  approximation  and 
possibly  appealing  to,  for  example,  the  method  of  multiple  scales  to  give  the  approximation 
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some  justification,  see  for  example  Floryan  and  Saric  (1979).  However,  it  is  now  accepted 
that  no  such  justification  can  be  made  and  that,  if  one  is  interested  in  the  evolution  of 
vortices  with  0(1)  wavelength  in  a  boundary  layer,  the  governing  equations  are  indeed  the 
parabohc  equations  solved  by  Hall  (1983),  thus  no  reduction  to  a  set  of  ordinary  differential 
equations  is  possible. 

In  the  nonlinear  regime  the  nonparallel  characteristics  of  Gortler  vortices  are  main¬ 
tained  and  the  disturbance  equations  must  be  solved  numerically.  Hall  (1988).  However 
in  the  small  wavelength  nonlinear  regime  much  analytical  progress  can  be  made  and  it 
is  possible  to  describe  vortices  so  large  that  they  have  an  0(1)  effect  on  the  mean  flow. 
Indeed  in  that  regime  it  is  found  that  the  mean  flow  is  actually  driven  by  the  vortices  over 
most  of  the  flow  and  adjusts  itself  so  as  to  make  all  small  wavelength  vortices  neutrally 
stable,  see  Hall  and  Lakin  (1988).  At  even  larger  amplitudes  the  strongly  nonlinear  states 
described  by  Hall  and  Lakin  (1988)  become  unstable  to  wavy  vortex  modes  as  discussed 
by  Hall  and  Seddougui  (1989). 

The  above  description  of  Gortler  vortex  growth  in  the  linear  and  nonlinear  regimes 
also  apphes  to  compressible  flows  it  the  Mach  number  is  not  too  large.  Hall  and  Malik 
(1989),  Wadey  (1990).  In  the  hypersonic  limit  Hall  and  Fu  (1989a,b)  showed  that  the 
linear  development  of  Gortler  vortices  becomes  much  simpler.  In  particular,  for  either  a 
Sutherland  or  Chapman  law  fluid,  the  first  term  in  the  expansion  of  the  neutral  Gortler 
number  in  terms  of  Mach  number  and  wavenumber  is  independent  of  nonparallel  effects. 
That  property  remains  true  at  higher  order  for  a  Chapman  Law  fluid  but,  for  a  Sutherland 
law  fluid,  nonparallel  effects  dominate  at  second  order. 

The  results  discussed  above  are  broadly  in  agreement  with  experimental  observations, 
certainly  the  nonparallel  incompressible  calculations  are  much  more  in  line  with  experi¬ 
mental  observations  than  are  the  parallel  flow  calculations.  Nevertheless  the  nonparallel 
theory  causes  philosophical  problems  for  the  transition  prediction  industry,  the  reason  why 
this  should  be  the  case  is  obvious.  Thus  in,  for  example,  the  e"  rule  for  transition  predic¬ 
tion  it  is  necessary  to  know  the  unique  growth  rate  for  a  vortex  at  a  given  downstream 
position,  but  the  main  result  of  the  nonparallel  work  is  that  no  such  quantity  exists.  In 
the  absence  of  a  unique  growth  rate  the  method  cannot  be  used  so,  in  some  cases,  the 
outcome  of  the  nonparallel  work  has  been  ignored  and  parallel  flow  theory  used  to  predict 
the  required  unique  growth  rate.  In  fact  an  unstable  Gortler  vortex  undergoes  most  of  its 
linear  growth  at  high  wavenumbers  so  the  latter  type  of  calculation  is  not  totally  flawed. 
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The  point  is,  of  course,  that  if  the  concept  of  a  unique  growth  rate  is  tenable  only  at  high 
wavenumbers,  and  as  the  growth  rate  can  be  found  in  a  simple  asymptotic  manner  there, 
the  high  wavenumber  theory  should  be  used  to  predict  the  growth  rate  there. 

The  motivation  for  much  of  the  research  in  recent  years  on  Gortler  vortices  has  come 
from  the  Laminar  Flow  Control  program  at  NASA  Langley.  In  particular  a  type  of  wing 
cross-section  developed  there,  see  Harvey  and  Pride  (1982),  has  significant  regions  of  cur¬ 
vature  on  its  underside.  These  concave  sections  are  required  to  control  flow  separation 
but  are  possible  causes  of  transition  via  the  Gortler  vortex  mechanism.  A  question  of 
some  importance  in  the  flow  around  this  type  of  wing  is  to  determine,  given  that  curva¬ 
ture  is  required  because  of  other  considerations,  how  the  curvature  should  be  distributed 
in  order  to  minimise  the  likelihood  of  transition  induced  by  Gortler  vortices.  A  matter 
of  equal  importance  is  that  of  whether  changes  to  the  curvature  distribution  made  to 
suppress  Gortler  vortices  wiU  enhance  the  likelihood  of  transition  being  caused  by  other 
mechanisms.  Further  motivation  to  study  compressible  Gortler  vortices  comes  from  the 
necessity  for  engineers  to  understand  the  flow  around  turbine  blades  or  the  flow  in  engine 
inlets.  In  fact  the  stability  problem  associated  with  the  latter  flow  is  made  much  more 
difficult  because  of  shock  waves  present  in  the  flowfield,  as  yet  no  progress  has  been  made 
with  the  Gortler  vortex  problem  in  the  presence  of  shocks  though  some  progress  has  been 
made  with  understanding  the  effect  of  shocks  on  travelling  wave  instabilities,  Cowley  and 
Hall  (1990). 

In  the  practical  situations  where  the  Gortler  mechanism  is  thought  to  be  important  it 
remains  an  open  question  if  free-stream  disturbances  or  wall  roughness  will  cause  the  initial 
vortex  growth.  Without  doubt  in  both  the  wind-tunnel  and  flight  situations  we  can  expect 
that  streamwise  vorticity  impinging  on  the  leading  edge  or  generated  by  imperfections  on 
the  wing  surface  will  be  present.  The  main  aim  of  the  present  paper  is  to  determine  the 
efficiency  of  surface  imperfections  in  generating  vortices.  This  wiU  be  done  over  the  whole 
range  of  Gortler  numbers,  this  means  that  we  will  have  to  discuss  the  spatial  inviscid 
instability  problem  for  Gortler  vortices,  apparently  this  has  not  been  done  previously. 

We  shall  see  that  the  investigations  of  Hall  (1982a)  at  high  Gortler  numbers  over¬ 
looked  a  wavenumber  regime  where  the  fastest  growing  spatial  Gortler  vortices  occurs.  A 
significant  property  of  these  modes  is  that  they  are  wall  modes,  actually  they  are  localized 
in  an  asymptotically  thin  layer  at  the  wall.  This  distinguishes  them  from  the  neutral  form 
of  the  small  wavelength  modes  of  Hall  (1982a)  which  are  trapped  in  a  thin  asymptotically 
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small  layer  in  the  interior  of  the  flow.  The  new  modes  which  we  describe  are  viscous  but 
they  connect  in  a  simple  way  to  the  in  viscid  spatial  modes  in  one  limit,  and  with  the  neu¬ 
tral  modes  of  Hall  (1982a)  in  another  limit.  Because  these  most  unstable  modes  choose  to 
locate  themselves  at  the  wall,  it  turns  out  that  they  couple  strongly  with  any  forcing  due 
to  a  surface  imperfection.  Our  calculations  suggest  that  the  stimulation  of  these  modes 
will  be  the  most  likely  consequence  of  a  quite  general  surface  imperfection,  therefore  they 
are  likely  to  be  the  source  of  transition  in  curved  flows  at  high  Gortler  numbers.  However, 
in  addition  to  these  modes,  we  shall  discuss  all  other  possible  types  of  induced  longitudi¬ 
nal  vortex  scalings.  Our  discussion  applies  to  walls  with  small  amplitude  imperfections  of 
somewhat  arbitrary  character  but  we  concentrate  mostly  on  the  case  when  the  perturba¬ 
tion  is  localized  in  the  spanwise  direction.  For  such  walls  we  find  that  the  vortices  grow  in 
a  wedge  shaped  region  downstream  of  the  obstacle,  though  in  the  immediate  downstream 
vicinity  of  the  obstacle  the  disturbance  decays. 

The  procedure  adopted  in  the  rest  of  this  paper  is  as  follows:  in  §2  we  formulate 
the  receptivity  problem  for  Gortler  vortices  induced  by  surface  imperfections.  In  §3  we 
consider  the  limiting  form  of  this  receptivity  problem  appropriate  to  the  small  wavelength 
limit;  in  this  limit  we  are  able  to  derive  an  asymptotic  form  for  the  coupling  co-efficient 
between  the  forcing  and  the  induced  vortex  field.  In  §4  the  0(1)  wavelength  problem  is 
investigated;  here  we  allow  the  forcing  to  vary  on  a  lengthscale  comparable  to,  and  shorter 
than,  the  boundary  layer  scale  in  the  flow  direction.  In  §5  we  isolate  the  fastest  growing 
mode  of  the  Gortler  stability  equations.  In  §6  the  receptivity  problem  for  this  mode  is 
discussed;  finally  in  §7  we  draw  some  conclusions. 

2.  Formulation  of  the  forced  Gortler  vortex  problem. 

We  consider  the  flow  of  a  viscous  fluid  over  a  wall  of  variable  curvature.  We  assume 
that  L  is  a  typical  lengthscale  over  which  the  curvature  changes.  If  is  a  typical  flow 
velocity  a  great  distance  from  the  wall,  and  v  is  the  kinematic  viscosity,  we  define  a 
Reynolds  number.  Re,  by 

(2.1) 

u 

and  we  suppose  that  with  respect  to  Cartesian  axes  x*  ,y*  ,z’  the  wall  is  defined  by 

y*  ^  L{Re-^g{x*/L)  +  ARe-^f{x*/L,Re^z*/L)}.  (2.2) 
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where  A  is  a  small  constant. 
We  define  variables  {x,y,z)  by 


{x*,Re3y*,Reiz*) 

L 


and  a  corresponding  velocity  vector  by 


.  4.  4.  4..  {u*,Re»v*,Re»w*) 

^00 

We  restrict  our  attention  to  the  limit  Re  00  and  write 

p"*"  =  p{x)  +  ARe~  » p{x,  y,  z)  +  0(  A^) 
(u+,i>+,u;+)  =  (u,V,0)  +  A(u,v,w)  +  O(A^) 


(2.3) 


(2.4) 


(2.5a, 6) 


where  is  the  pressure,  scaled  on  pU^ ,  with  p  the  fluid  density.  Moreover,  11  and  V 
depend  only  on  x  and  y  whilst  u,v,w  depend  on  all  three  dimensionless  co-ordinates.  We 
first  take  the  limit  A  — »  0  with  Re  held  fixed  and  equate  terms  of  order  A®,  A  in  the 
Navier  Stokes  equations.  In  the  resulting  systems  for  (Q,O,0,p)  and  (u,v,u),p)  we  take 
the  further  limit  iZe  — ♦  00  to  obtain,  at  the  zeroth  order  level  of  approximation  in  Re, 


G.  +  Dv  =  0, 

GH,  -t-  GGy  =  — p,  -t-  Gyy, 

and 

u,  -f  u,  -I-  =  0. 

UU,  +  OUy  -I-  UG,  -|-  VUy  =  A2U, 

Gu,  +  VVy  -t-  GG,  -I-  VVy  =  -Py  -f-  A2V, 
UW,  +  VWy  =  — P,  -1-  A2U>- 

Here  the  operator  A2  is  defined  by 


(2.6) 


(2.7) 


^2  = 


and  by  writing  down  a  Taylor  series  expansion  for  the  no-slip  condition  about  y  =  g  we 
obtain  the  conditions 


G  -  G  =  0,  y  =  g{x), 

u=-fUy,  v  =  -fVy,  w  =  0,  y  =  g{x). 


(2.8a, 6) 
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Meanwhile  for  large  values  of  y  we  require  that 


U  — >  u,(z),  y  —*  oo 
{u,v,w)-^0,  y -*  oo, 


(2.9a, 6) 


where  u,{i)  is  the  dimensionless  free-stream  velocity;  in  this  paper  we  shall  confine  our 
attention  to  Blasius  flow  so  that  u,  =  1.  In  order  to  bring  out  more  clearly  the  destabilizing 
effect  of  wall  curvature  we  make  the  Prandtl  transformation 


y-^y-9, 

c  -♦  0  +  g'n, 

5  — >  t>  4-  g'u. 

whilst  all  other  variables  remain  unchanged.  In  this  case  (2.6)  -  (2.9)  become 


fix  +  Cy  =  0, 

QQ,  +  VUy  =  — P,  4-  Uyy, 

0  =  0  =  0,  y  =  0, 

Q  — >  u,(z),  y  — »  oo. 


(2.10) 


and 


U,  Vy  +  Wt  =  0, 

Qu,  4-  ^Uy  4-  vUy  4-  uu,  =  A2U, 

QC,  4-  VVy  4-  uO,  4-  vVy  4-  Gx{x)nu  =  -py  +  A2V, 

(2.11) 

UWg  4-  ^Wy  =  — p,  4-  A2«>, 

u  =  —fuy ,  i)  =  0,  ti)  =  0  y  =  0, 


u,v,  w  —*  0,  y 


00. 


drg 

In  (2.11)  we  have  replaced  2^-2  by  Gx{x)  where  G,x  will  be  referred  to  as  the  Gortler 
number  and  wall  curvature  respectively.  Thus  (2.10)  must  first  be  solved  for  the  basic 
state  whilst  (2.11)  then  determines  the  forced  vortex  flow.  We  further  note  that  the  wall 
condition  on  v  could  be  modified  so  as  to  determine,  for  example,  the  effect  of  wall  suction 
on  the  generation  of  vortices.  Finally  in  this  section  we  assume  that  /  may  be  written  as 


/  = 


6 


so  that  (2.11)  can  be  Fourier  transformed  in  z  to  give 

u,  +  V,  +  ikw  =  0, 

2 

GU,  +  CUy  +  UGy  +  UG,  =  Uyy  —  fc  U, 

QVg  +  VVjf  +  VVy  +  UVg  +  GXGU  =  — Pv  +  Vyy  — 

(2.12a  -  /) 

Giy,  +  dWy  =  — tfcp  +  Wyy  —  H  Wy 

u  =  q{k)F{x)y  V  =  0,  to  =  0,  y  =  0, 
u,v,w  — ♦  0  y  — »  oo. 

Here  k  is  the  transform  variable,  —q  is  the  transform  of  q,  u  is  the  transform  of  u,  etc.,  and 
we  have  defined  F{x)  =  Gy(x,0)/(x).  In  the  following  sections  we  discuss  the  solutions  of 
(2.12)  in  a  variety  of  situations;  for  convenience  we  note  that  if  p  and  w  are  eliminated 
from  the  disturbance  equations  we  obtain 

{  U,yy  +k^+  k^Vy  }  U  +  G,  Uy  y  +  {  G , ,  y  +  U 

d 

+  +  2{G.y  +  C,— }u,  (2.12y) 

+  Vyyyy  “  GUyyy  "  {  Gy  +  2fc^}Uyy  +  {G,y  +  fc  ^  G }  Gy  =  0. 

Thus  the  problem  for  the  forced  vortex  flow  can  be  reduced  to  the  solution  of  (2.12a,g) 
subject  to  the  boundary  conditions  (2.12e,f).  At  finite  values  of  G  the  solution  of  this 
system  can  only  be  found  by  numerical  integration.  However  it  is  first  instructive  to 
consider  the  receptivity  problem  for  vortices  with  a  wavelength  appropriate  to  the  right 
hand  branch  of  the  neutral  curve;  the  results  we  obtain  are  invaluable  in  interpreting  the 
calculations  at  0(1)  wavelength  which  will  be  discussed  in  §4. 

3.  The  small  wavelength  limit. 

We  now  consider  the  asymptotic  limit  fc  — >  oo  in  which  case  the  unforced  problem 
becomes  unstable  in  a  layer  of  thickness  k~  *  located  in  a  position  which  maximises  the 
downstream  growth  of  the  vortex.  From  Hall  (1982a, b)  we  know  that  small  wavelength 
vortices  develop  in  a  nonparallel  manner  in  a  fc"  ^  neighbourhood  of  the  neutral  location. 
Thus,  if  the  Gortler  number  G  is  written  as 

G  =  (3.1) 
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then  suppose  that,  using  the  approach  of  Hall  (1982a),  we  find  that  correct  to  zeroth  order 
in  k~  ^  the  neutral  location  is  given  by  a:  =  I.  We  then  write 

X  =  k{x-  X).  (3.2) 

We  will  allow  the  wall  forcing  function  F{x)  to  be  slightly  more  general  than  that  assumed 
in  §2;  more  precisely,  we  now  allow  F  to  vary  on  the  X  lengthscale  so  that 

F{x)  =  TiX).  (3.3) 

We  will  now  show  how  the  wall  roughness  induces  a  longitudinal  vortex  structure  within 
a,  k~^  neighbourhood  of  x  =  Z;  notice  that  if  we  wish  to  recover  the  case  with  F  varying 
on  X  lengthscale  we  simply  need  to  replace  F'{X)  by  a  constant  in  the  following  analysis. 

It  follows  from  (2.12)  that  viscous  effects  in  the  normal  and  spanwise  directions  are 
comparable  when  ^  ~  0(fc);  this  is  achieved  in  the  bulk  of  the  flow  using  a  WKB  type  of 
expansion.  However,  as  discussed  by  HaU  (1982a),  this  expansion  fails  near  the  wall  and 
where  UUy  has  a  local  maximum.  We  suppose  that  this  maximum  occurs  at  y  =  ^  so  that 
the  forced  flow  must  have  a  region  II  of  depth  k~  >  centred  on  y  =  y.  Above  and  below 
this  region  are  two  further  regions,  I  and  III,  where  the  WKB  structure  is  appropriate. 
At  the  bottom  of  region  III  the  WKB  structure  will  be  seen  to  break  down;  in  fact  the 
six  different  WKB  exponents  collapse  into  just  two  distinct  values  as  y  — >  0.  Thus  an 
adjustment  or  transition  layer  is  needed  as  y  — »  0;  we  shall  see  below  that  this  is  of  depth 
k~*  and  we  refer  to  this  region  as  IV.  Finally,  where  y  =  0(fc“^),  the  viscous  derivatives 
in  the  y  direction  again  become  important  so  that  a  wall  layer,  V,  of  depth  k~^  is  needed. 
The  different  regions  are  illustrated  in  Figure  (3.1). 


The  wall  layer  solution 

Here  we  define  ^  by 

^  =  ky 

and  the  wall  conditions  suggest  that  u  and  v  are  0(1)  here.  In  fact  it  can  be  seen  from 
(2.12a)  that  if  u  is  0(1)  in  V  then  v  must  be  of  order  k;  we  are  therefore  led  to  the  expansion 


u  =  uo(0  +  •••. 

V  =  kvo{^)  +  •••■, 


(3.4a,  b) 
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and  we  expand  the  basic  velocity  component  G  in  this  layer  as 


C  =  A^A:-^  + 


(3.5) 


where  X{x)  is  the  wall  shear  at  x  =  2.  A  similar  expansion  can  be  written  down  for  V 
but  is  not  needed  here.  If  we  substitute  the  above  expansions  into  (2.12)  we  find  that  the 
zeroth  order  approximation  to  that  system  yields 


2 

{ :772  -  ^0  = 

no  =  ^,VQ  =  vof  =  0,  ^  =  0. 

The  solution  of  (3.6)  which  does  not  grow  exponentially  with  ^  is  then  seen  to  be 

no  =  J^{X)e-i,VQ  =  +  3^2}. 


so  that  when  ^  »  oo  (3.7)  gives 

u  ~  e  + 


V  ^ 


-XGXJ^k^ 

24 


iyk*) 


i\3._ 


*v 


+ 


(3.6) 


(3.7) 


(3.8) 


Solution  in  IV 

We  can  see  from  (3.8)  that  the  next  region,  IV,  is  of  depth  k~  *  so  we  put 


(3.9) 


and  the  appropriate  expansions  of  u  and  v  are  then 


u  =  UQ{Ti)e  **  +  ••• 
v  =  kiVo{T))e-^^  +  ■■■ 


and  after  substituting  into  (2.12) 
equations  we  find 

d^Up 

drj^ 

um 


and  eliminating  Vp{r])  from  the  resulting  zeroth  order 


xGX'^ri 
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Up  = 


0 


=  C/6'  =  0,  U{Q)  =  T{x). 


(3.10) 
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The  solution  of  this  system  is 


where  a  =  ‘  and  satisfies 


Uq  = 


(3.11) 


-^-c^  =  0,  «'(0)  =  1,  ^''(0)  = ’J'"(0)  =  0. 

For  large  values  of  C  the  three  possible  forms  for  ^  have 

~  e*^  ,  e«*  ^ 

respectively  so  that  grows  exponentially  when  C  —*  oo.  In  fact  the  asymptotic  form  of 
^  for  large  C  is  found  to  be 

/3 

~  ^exp  \^-C*  j  +  ••• 

where  w  is  a  constant  given  to  three  decimal  places  by  u;  =  0.285. 

It  follows  that  for  large  values  of  r?  we  can  write  u  as 


(3.12) 


In  fact  the  two  other  (complex)  exponentially  decaying  solutions  match  onto  exponentially 
decaying  WKB  solutions  in  III  which  do  not  have  a  turning  point  layer  at  j/  =  g.  Thus  it 
is  sufficient  for  us  to  consider  the  matching  of  (3.12)  with  the  WKB  solution  which  does 
indeed  have  a  turning  point  layer  at  y  =  p;  it  is,  of  course,  the  part  of  the  velocity  field 
which  drives  the  potentially  unstable  vortex  velocity  field. 


The  WKB  solution  in  III 

In  this  region  we  seek  asymptotic  solutions  for  u  and  v  which  take  the  form 

(u,u)  =  e‘'^o  +  •  •  •}.  (3.13) 

Here  9{y)  is  the  WKB  phase  function  and  the  relative  scalings  of  u  and  v  follow  from  the 
requirement  that  {d^  +  ~  FOy  in  III;  the  latter  balance  is  the  crucial  one  shown 

by  Hall  (1982a)  to  be  necessary  for  small  wavelength  vortices.  It  is  a  routine  matter  to 
substitute  for  {u,v)  from  (3.13)  into  (2.12)  and  solve  for  (uoi^’o)-  In  fact  at  zeroth  order 
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uo  and  DO  are  not  determined  but  the  consistency  of  a  pair  of  linear  equations  for  these 
functions  yields  the  phase  equation 


(02_  1)3  ^ 


where  u[y)  —  u{x,y).  In  order  that  u  and  v  should  decay  with  y  we  must  reject  values  of 
9  with  positive  real  part.  The  three  acceptable  solutions  are  then  found  to  be 


0  = 

^  =  -  [1  - 


(3.14a, 6) 


where  the  branch  of  the  square  root  in  the  right  hand  half  plane  is  to  be  taken.  The  first 
of  these  roots  vanishes  where  1  =  Gx^Uy  and  this  necessitates  the  existence  of  the  turning 
point  layer  II.  The  complex  roots  lead  to  no  such  turning  point  layer  so  that  the  part  of 
the  WK3  solution  appropriate  to  these  roots  is  valid  in  0  <  y  <  oo.  In  fact  these  solutions 
match  directly  onto  the  part  of  the  solution  in  III  appropriate  to  the  complex  valued  large 
solutions  of  (3.12).  Thus  we  neglect  this  part  of  the  WKB  solution  since  it  does  not 
interact  in  any  way  with  the  induced  vortex  instability  which  we  will  calculate  in  II.  We 
therefore  write  u  in  III  as 


_-k  /*  [1- {Gxuuv}  1)3  dv 


{uQ  +•••}  +  E.D. 


where  E.D.  denotes  the  other  WKB  solutions  not  having  a  turning  point  layer  at  y  =  y. 
At  higher  order  the  amplitude  function  no(A',y)  is  found  to  satisfy 


—  [0^  —  iJGA"  G  (y)uo  —  2GuUyXv0uo  —  2Gx^^^v^0x 

+  Q9Gx^%UQy  +  30'(0^  —  1)^(1  —  50^)uo  —  AOxGuUyyUQ  =  0. 


Here 

e  (y)  =  [wuyxlxj 

I  X  =z 

and  V  =  d{x,y).  Thus  no  satisfies  a  partial  differential  equation  with  coefficients  dependent 
on  both  X  and  y.  Non-parallel  effects  lead  directly  to  the  term  proportional  to  X  in  (3.15). 
Furthermore  in  order  to  match  with  the  solution  in  IV  we  require 


no 


k  *  [xGA^]  13 


■E(x),  y  ->  0. 


(3.1G) 
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Ill  order  to  siiiiplily  the  solution  of  the  equations  for  uq  we  w’lite 


where 


with 


and 


UO 


2Tu;(5Jt[/o 


—  1^1  ^exp{X]  {x,y)}, 


yl'T[xGA2]tV 


7  =  rH4?#} 


N{y)  =  ixuUy)  -  l)(xuUy)*j  +  4u^Uj,X7} 


»-^iy  — y 


=  2  = 


(3.17) 

(3.18) 

(3.19) 

(3.20) 


We  note  at  this  stage  that  N  {y)  =  0  so  that  the  integral  /q  ^  dy  converges  when  y  =  y, 

even  though  d{y)  =  0.  After  writing  uq  in  the  form  (3.17)  we  find  that  Uo{x,y)  satisfies 


OUq  u  OUq 

dy  39  dX  ’  (3.21) 

Uo  =  y  =  0. 


This  system  is  most  easily  solved  by  taking  the  Laplace  transform  of  Uq.  However  before 
taking  this  transform  it  is  convenient  to  suppose  that  ^(X)  =  0  for  X  <  X'  and  that 
u  =  u  =  0  for  AT  <  X' .  We  now  make  the  change  of  variable  X  —  X  —  X'  and  let  UQ{s,  y) 
be  the  transform  of  Uq.  We  then  find  that 

Uq  =  f{sy -^0  (3.22) 


where 

f{s)=[  +  A'*)dX.  (3.23) 

We  must  now  find  the  limiting  form  of  the  WKD  solution  when  y  —*  y',  this  will  provide 
us  with  the  required  matching  condition  which  we  need  in  order  to  find  the  forced  motion 
in  II. 

Firstly  we  note  that  when  y  y  the  phase  function  9  may  be  written  as 


where 


-*/)  +  •••  (3.24) 

(3,25) 
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and  a  *  denotes  a  function  evaluated  at  a:  =  f,t/  =  p.  It  follows  that  in  the  limit  y  y 


7^  ~  -^\y  -  y\^-  ^  expiXl  +  +  k  ['  Ody  + 

1 0  A;  T 


k{y  -  g)'^ 


M},  (3.26) 


with 


^  _  2«  (u*)j  exp{jj/Q*'  vdy  +  /p  jjl/p  ^^)dy} 
MT[xGA2ltT. 


(3.27) 

Here  /  denotes  the  finite  part  of  an  integral.  Finally  we  note  that  Uq  in  the  limit  y  g 


takes  the  form 

Uq  ~  f{^)\y  -  g\^e‘^, 

where 


(3.28) 

(3.29) 


Solution  in  the  transition  layer 

Here  the  solution  is  essentially  that  found  by  Hall  ( 1982a, b)  so  that  we  only  summarize 
briefly  how  it  is  obtained  from  the  disturbance  equations.  The  layer  is  of  depth  k~  a  so  we 
define  Y  by 

Y  =  k^y-g} 

and  then  u  and  o  e.xpand  as 

u  =  Uo{Y)  +  k-^Ui{Y)  +  ■■■ 
k^Vo{Y)  +  kH\{Y)  Y  ■■■. 

At  zeroth  order  we  obtain  the  equations 

Uq  +  HoS*  =  0, 

+  \Gu'Uo  -  0, 

which  have  a  consistent  solution  if 

xGa*  a*  =  2. 

Thus  at  zc'ioth  order  the  leading  order  disturbance  velocity  field  remains  unknown.  In  fact 
since  y  was  chosen  such  that  (u^)yy  vanishes  at  y  =  y’  we  must  proceed  to  the  0(A:”^) 
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appioxiiiiutioii  level  to  the  disturbance  equations  before  Uq  and  Vq  are  found.  At  that 
order  we  find  that  Uq  satisfies 


O'^Uq 

OY^ 


20*  dUo 
~~dX 


+  -jixGuu^yiUo  - 


A/2y2 

4 


Uo  =  0 


(3.30) 


In  the  absence  of  any  forcing  we  would  seek  solutions  of  this  evolution  equation  which  have 
Uq  — >  0,  |y|  — >  oo.  However,  for  the  forced  problem  we  replace  the  condition  at  Y  =  — oo 
by  one  which  enables  us  to  match  with  the  solution  in  III.  With  that  in  mind  we  take  out 
the  X  dependence  of  Uq  appropriate  to  (3.26)  by  writing 


Uo  =  Uo{X,Y)e^^-^^" 


in  which  case  Uq  satisfies 


O^Uq  2_,dUQ 
OY'^  ~  3^  dX 


2v'2 


M^Y 


2=. 


Uq  —  -u  IUq  =  0. 
4  3 


We  can  then  make  the  change  of  variable  X  =  X  —  X'  and  take  the  Laplace  transform 
of  this  equation  with  respect  to  X .  If  we  denote  the  Laplace  transform  of  Uq  by  Uq  we 
find  that  Uq  satisfies  a  parabolic  cylinder  equation.  More  precisely,  if  Uq  is  to  vanish  when 
Y  — »  oo,  we  require  that 

Uq  =  BU{'^{s  +  I),MrY)  (3.31) 

where  U{a,  x)  is  a  parabolic  cylinder  function.  At  this  stage  if  s  is  chosen  such  that 


2(2* 

3M 


[s  +  I] 


—n  — 


1 

2 


then  Uo(±oo)  =  0.  This  choice  of  s  corresponds  to  the  linear  eigensolutions  found  by 
Hall  (1982a, b).  In  order  to  match  the  forced  solution  with  that  in  III  we  note  that  when 
Y  — >  — oo 


Uq 


D\/^ 


+  -^1  +  2) 


Thus  matching  with  the  coreflow  is  achieved  by  choosing  B  =  B{s)  such  that 


+  -^1  +  2) 


f +  3 ^ +  fo^ [-^0  ^^y^ ^y) 


Having  determined  D  we  can  write  down  the  inverse  transform  for  uq  and  the  flow  in 
II  is  then  completely  known.  However,  the  gamma  function  has  a  sequence  of  poles  at 
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the  origin  and  on  the  negative  real  axis  so  that,  if  /  has  no  poles  or  branch  cuts  Uq  will 
be  found  from  the  contribution  of  these  poles.  The  least  stable  mode  corresponds  to  the 
pole  at  the  origin  so  that  for  large  X  the  solution  will  be  dominated  by  that  mode.  As  a 
measure  of  the  size  of  the  disturbance  in  the  layer  III  we  can  take  the  centreline  velocity. 
After  some  manipulation  we  see  that  the  streamwise  disturbance  at  V'  =  0  for  large  values 
of  X  is  given  by 


4.2“  •  O' 


u  ~ 


exp{y{X*  +  Xf  +  IX*  - -^X} 


where 


It  follows  from  the  above  discussion  that  the  coupling  coefficient  which  relates  a  typ¬ 
ical  vortex  velocity  field  with  the  forcing  function  is  proportional  to  k*  exp{— fc/g  |1  — 
{Gx^Sv}’ Thus,  since  this  result  is  valid  only  in  the  limit  fc  1  we  see  that  small 
wavelength  wall  roughness  is  an  extremely  inefficient  generator  of  Gortler  vortices.  The 
reason  why  this  should  be  so  is  evident  from  the  above  discussion;  it  is  of  course  a  direct 
consequence  of  the  fact  that  the  vorticity  induced  by  the  roughness  must  diffuse  over  an 
0(1)  distance  to  the  layer  where  potentially  unstable  vortices  of  this  wavelength  can  ex¬ 
ist.  We  conclude  that  small  wavelength  vortices  of  wavenumber  proportional  to  G«  are 
almost  certainly  unlikely  to  be  generated  by  wall  roughness.  Thus  it  would  be  extremely 
surprizing  if  small  wavelength  vortices  could  be  induced  experimentally  by  wall  roughness; 
however  the  work  of  Hall  (1990)  shows  that  such  vortices  could  be  induced  by  free  stream 
disturbances.  We  further  note  that,  whatever  the  initial  size  of  the  vortex  wavelength,  the 
fact  that  it  is  observed  experimentally  that  the  vortex  wavelength  is  conserved  as  it  devel¬ 
ops  in  the  x  direction  means  that  ultimately  the  small  wavelength  limit  becomes  relevant. 
Now  we  shall  relax  the  constrain  on  the  size  of  the  forced  vortex  wavelength  and  see  if  a 
stronger  coupling  between  the  induced  form  and  the  forcing  is  possible. 


4.  Vortices  of  0(1)  wavelength 

Here  we  will  be  primarily  interested  in  the  case  when  the  wavelength  of  the  forced 
vortex  is  comparable  with  the  boundary  layer  thickness.  We  shall  concern  ourselves  with 
both  ‘isolated’  and  ‘distributed’  roughness  elements.  For  convenience  we  shall  initially  take 
q  =  I  which  corresponds  to  a  delta  function  shaped  hump  in  the  spanwise  direction,  later 
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we  choose  q  to  correspond  to  other  shapes.  Moreover  we  shall  also  discuss  the  situation 
when  the  wall  forcing  function  varies  on  a  fast  lengthscale,  for  an  isolated  roughness  element 
this  type  of  forcing  provides  unique  initial  conditions  for  the  partial  differential  equations 
describing  the  vortex.  Later  in  this  section  we  shall  return  to  the  case  when  F{x)  varies 
on  an  0(1)  lengthscale,  we  shall  see  that  this  situation  is  the  most  efficient  in  producing 
linearly  growing  Gortler  vortices  of  0(1)  wavelength. 

We  shall  first  determine  solutions  of  (2.12)  appropriate  to  the  case  when  F{X)  varies 
on  a  relatively  fast,  0(f),  lengthscale.  We  suppose  that  the  forcing  begins  at  x  =  2  and 
write 


e 


and  then 

F{x)  =  F*iX).  (4.1) 

In  fact  without  any  loss  of  generality  we  can  take  J  j  so  ;hat  the  original  lengthscale  L 
has  then  been  fixed  in  terms  of  the  distance  from  the  leading  edge  to  the  position  where 
the  forcing  begins.  Clearly  we  expect  a  different  response  of  the  flow  when  S  is  before  or 
after  the  unstable  regime  of  a  vortex  with  wavenumber  k\  with  J  now  fixed  we  must,  of 
course,  investigate  that  possibility  by  varying  the  Gortler  number  G. 

In  order  to  find  the  forced  flow  in  a  neighbourhood  of  £  we  first  note  that  when 
y  — >  0,  nd,  ~  y/e  so  that  convective  and  vertical  diffusion  effects  are  comparable  in  a  layer 
of  depth  €» ,  and  the  wall  forcing  implies  that  u  will  be  0(1)  there.  It  follows  that,  if  u  is 
indeed  0(e°)  in  this  layer,  then  v  will  be  0(f'»)  there.  We  deduce  from  (2.12c,d)  that  p 
and  w  must  then  be  0(e“  » ),  0(e“  ^)  respectively.  Thus  if  we  put 

C  =  y/f* 

then  for  ^  =  0(1)  we  seek  a  solution  of  (2.12)  which  takes  the  form 

{u,v,w,p)  =  {uo{X,^),  e"«i;o(-^,0>  *  Po{X  ,^))  +  ■  ■  ■ 


whilst 


where  A  =  Uy{^,0),p 
find  that 


(!2,0)  =  +  ■■■  (4.20,6) 

2Vyy{^,0).  If  we  substitute  the  above  expansion  into  (2.12g) 


we 


=  0, 


(4.3) 
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where  5  is  the  Laplace  transform  variable  and  Do  is  the  transform  of  uq.  If  t)0{  is  to  be 
finite  at  infinity  then 

^=aJ  MX>sU)d^  +  B.  (4.4) 

a?  •'O 

Here  A  and  B  are  constants;  however,  the  transformed  continuity  equation  evaluated  at 
the  wall  yields 

B  =  -sF{s).  (4.5) 

Here  F  is  the  transform  of  F*  and  then  vanishes  at  infinity  if 

A=—3s*\*B.  (4.6) 


(4.11) 


Thus  Uq  and  Fq  decay  to  zero  exponentially  as  y  — >  oo  and  satisfy  the  matching  conditions 
(4.8).  In  each  layer  the  transforms  can  be  inverted  and  the  large  X  form  of  this  velocity 
field  can  be  used  as  an  initial  condition  for  the  full  linearized  disturbance  equations.  Now 
we  consider  the  special  case  F*  =  6{X)  where  6  is  the  delta  function.  It  is  easy  to  show 
by  inverting  (4.7)  and  the  corresponding  form  for  uq  that  this  choice  of  F*  leads  to  a 
similarity  solution 

UQ'^l/XuaiHX*),  VQ'^lIX^VQiilX*). 

In  fact  the  above  similarity  solution  is  also  the  large  X  asymptotic  form  of  a  disturbance 
velocity  field  induced  by  a  forcing  function  F*{X)  of  compact  support.  We  conclude  that 
the  similarity  solution  can  be  used  as  the  initial  condition  for  the  disturbance  equations 
in  the  case  when  vortices  are  induced  by  an  isolated  roughness  element.  We  note  here  a 
related  similarity  solution  situation  arises  in  linear  triple  deck  theory,  Hunt  (1971),  Smith 
(1973),  but  here  the  absence  of  a  streamwise  pressure  gradient  in  the  Gortler  equation 
causes  the  velocity  field  not  to  be  confined  to  the  wall  layer.  We  will  now  show  how  the 
similarity  solution  can  be  derived  directly  from  the  disturbance  equations. 


The  Similarity  Solution. 

The  appropriate  similarity  variable  implied  by  the  above  discussion  is 

C  =  y(^)*  (4.12) 

X 

where 

X  =  I  —  I  (4.13) 

and  for  C  =  0(1)  the  disturbance  velocity  components  u  and  v  are  found  to  expand  as 


At  _ 

u  =  — wo(C)  +  •  •• 

I  (4.14a, 6) 

V  =  — vo(C)  +  ••• 

X* 

and  at  the  wall  we  now  impose  the  conditions 


uo  =  u'o  =  uo  =  0,  C  =  0- 


(4.15) 
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The  functions  uo  and  vq  satisfy 


where  q*  is  as  defined  by  (4.11).  In  Figure  (4.2)  we  have  shown  the  function  q*{y,k) 
for  several  values  of  the  wavenumber  k.  We  note  that  for  the  smaller  values  of  k  the 
function  q*  has  a  maximum  value  away  from  the  origin.  At  larger  values  of  k  the  function 
q*  ~  e”**  and  therefore  decreases  monotonically  as  y  increases  from  zero.  We  can  then 
combine  (4.10),  (4.14)  to  give  a  composite  velocity  field  at  some  small  value  of  x  which  can 
be  used  as  initial  conditions  for  the  full  Gortler  disturbance  equations.  Before  discussing 
our  numerical  results  obtained  for  such  a  ceilculation  we  note  that  the  similarity  solution 
obtained  above  is  unique  only  up  to  a  multiplicative  constant.  That  constant  can  be 
evaluated  for  a  particular  wall  forcing  function  by  obtaining  the  large  X  form  of  the 
velocity  field  obtained  by  inverting  the  transformed  velocity  field.  The  constant  will  be 
0(1)  and  is,  of  course,  a  function  of  the  particular  function  chosen.  Moreover,  this  ‘coupling’ 
constant  gives  a  measure  of  the  input  disturbance  velocity  field  for  the  marching  problem. 
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We  see  below  that  at  the  position  where  the  downstream  velocity  field  induced  by  this 
disturbance  begins  to  grow  it  has  decreased  in  magnitude  several  orders  of  magnitude  so 
that  the  boundary  layer  is  extremely  unreceptive  to  this  type  of  forcing.  For  that  reason  it 
seems  unnecessary  for  us  to  gives  values  for  the  0(1)  coupling  constant  relating  the  forcing 
function  and  the  similarity  solution. 

Now  we  shall  report  on  some  results  we  have  obtained  by  integrating  the  Gortler  vortex 
equations  with  the  similarity  solution  discussed  above  used  to  provide  an  initial  vortex 
velocity  field.  A  composite  solution  was  used  to  begin  the  calculations  with  x  =  .5  +  Ax 
for  Ax  small.  The  results  presented  here  were  obtained  with  Ax  =  .001,  1600  points  in 
the  y  direction  and  a  vertical  step  length  =  .02.  The  results  were  of  course  verified  by 
varying  these  quantities  in  particular  cases.  It  is  of  interest  to  note  that  the  shape  of  the 
profile  q(y,  k)  for  small  k  shown  in  Figure  (4.2)  means  that  the  composite  normal  velocity 
field  will  have  a  double  maximum  because,  in  Figure  (4.1),  the  inner  velocity  field  also 
has  a  local  maximum.  We  further  note  that  for  fc  »  1,9*  ~  e"**.  The  scheme  used  to 
march  the  Gortler  equations  downstream  is  precisely  that  used  by  Hall  (1983)  so  we  give 
no  details  of  it  here.  As  a  measure  of  the  disturbance  energy  we  used  E{x)  defined  by 


and  the  local  growth  rate  (5  =  E~  Following  Hall  (1983)  the  position  of  neutral  stability 
can  be  defined  as  the  location  where  /3  =  0.  This  position  will  to  some  extent  depend  on 
the  choice  of  the  flow  property  used  to  define  0.  Our  calculations  for  several  values  of 
k,G  showed  that  E  decreases  by  several  orders  of  magnitude  before  P  becomes  positive. 
Effectively  this  means  that  the  type  of  localized  forcing  discussed  so  far  in  this  section  will 
generate  Gortler  vortices  but  the  vortices  produced  by  the  forcing  are  significantly  reduced 
in  size  before  they  begin  to  grow.  However,  in  the  absence  of  any  other  forcing,  this  type 
of  mechanism  will  be  that  responsible  for  the  growth  of  vortices  and  so  it  is  worth  giving 
some  more  detailed  results.  For  several  values  of  k  with  fixed  values  of  G  we  computed 
the  position  at  which  the  induced  vortex  structure  begins  to  grow.  The  local  wavenumber 
and  Gortler  number  were  calculated  at  that  point  and  used  to  generate  a  neutral  wave  in 
the  local  wavenumber  {kg),  local  Gortler  number  (G)  plane.  These  curves  are  plotted  in 
Figure  (4.3)  for  G  =  1,2,8.  It  is  important  to  realize  at  this  stage  that,  since  the  forcing 
was  imposed  at  the  same  value  of  x  in  each  case,  the  value  of  G  used  is  a  measure  of 
how  close  the  forcing  is  to  the  unstable  region.  In  fact  there  is  very  little  variation  in  the 
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different  curves  so  we  conclude  that  localized  forcing  generates  vortices  which  grow  at  a 
local  Gortler  number  of  about  10.  Here  it  is  of  interest  to  note  that  the  corresponding 
critical  Gortler  number  appropriate  to  free-stream  disturbances  is  about  6  so  that  free- 
stieam  disturbances  are  the  most  dangerous.  Indeed  the  energy  loss  of  a  vortex  induced 
by  a  free-stream  disturbance  is  only  an  0(1)  fraction  of  its  initial  value  so  it  would  appear 
that  the  localized  wall  roughness  mechanism  would  be  relevant  only  in  an  experimental 
facility  with  remarkably  small  free-stream  disturbances. 

In  Figures  (4.4a,b)  we  show  the  downstream  development  of  the  velocity  field  induced 
by  a  delta  function  disturbance  imposed  at  a:  =  .5  with  k  =  .45,  G  =  8.  We  notice  that 
the  initial  form  of  the  streamwise  disturbance  velocity  component  corresponds  to  a  double 
vortex  velocity  field  because  of  the  change  of  sign  of  u.  This  is  possibly  one  of  the  reasons 
why  localized  wall  roughness  is  not  an  efficient  generator  of  Gortler  vortices.  We  further 
note  that  the  velocity  field  at  the  larger  values  of  x  is  quite  typical  of  the  non-parallel 
velocity  fields  obtained  by  Hall  (1983).  In  order  to  quantify  the  inefficiency  of  localized 
wall  roughness  at  generating  Gortler  vortices  let  us  consider  Figures  (4.4b)  in  more  detail. 
The  local  boundary  layer  thickness  for  Blasius  flow  grows  like  x»  so  that  the  local  Gortler 
numbers  grow  like  xj.  Thus  when  x  increases  from  0.521  to  1.401  the  Gortler  number 
increases  from  about  8  to  40.  However,  even  in  this  highly  unstable  regime  the  maximum 
value  of  V  is  less  than  10%  of  its  initial  value.  Later  in  this  paper  we  shall  see  that  it  is 
possible  to  choose  the  wall  forcing  function,  and  the  vortex  wavelength,  to  produce  vortices 
with  an  almost  negligible  loss  in  amplitude. 

Now  we  shall  discuss  some  solutions  we  have  obtained  for  the  forced  vortex  problem 
when  the  forcing  occurs  on  an  0(1)  streamwise  lengthscale.  The  equations  (2.12b),  (2.12g) 
were  solved  using  the  approach  of  Hall  (1983)  subject  to  the  conditions  (2.12e,f).  Here  we 
shall  report  on  calculations  for  the  case  when  F{x),  the  forcing  function,  is  given  by 

F  =  40(x  - -)e-20(»-7)’,  x>i  (4.19) 

2  2 

The  disturbance  equations  for  a  given  k  were  marched  forward  from  x  =  0  with  zero  initial 
disturbance.  For  a  given  function  q{k)  we  can  invert  the  Fourier  transform  in  z  to  give 
the  velocity  field  induced  by  a  hump  of  height  proportional  to  q{z),  the  inverse  of  q{k). 
Before  discussing  our  results  we  point  out  that  further  calculations  for  wall  shapes  other 
than  (4.19)  were  carried  out  and  give  qualitatively  similar  results. 
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In  Figures  (4.5a,b)  we  show  the  development  of  the  functions  u  and  v  corresponding 
to  the  case  (4.19)  with  k  =  .45,  G  =  8.  We  note  that  at  smaller  values  of  x  the  maximum 
of  u  occurs  at  the  wall.  We  see  that  subsequently  the  maximum  value  of  u  occurs  at  a 
position  which  moves  away  from  the  wall  and  that  the  smallest  maximum  value  occurs  at 
about  X  =  3.5  and  is  approximately  .26.  This  is  to  be  compared  with  a  maximum  value 
of  4.84  of  F{x).  Thus  we  now  see  that  for  a  non  localized  wall  forcing  there  is  a  strong 
coupling  between  the  induced  vortex  field  and  the  forcing.  Furthermore  we  point  out  that 
the  ratio  between  the  smallest  maximum  of  u  and  the  maximum  of  F  is  a  function  of  k  and 
G.  Thus  by  ‘tuning’  k  and  G  we  can  in  principle  maximise  the  coupling  coefficient  between 
the  vortex  field  and  the  forcing.  However  such  a  calculation  was  not  carried  out  since  it 
would  require  a  great  deal  of  computer  time.  Moreover  the  results  of  Figures  (4.5a,b)  are 
sufficient  for  us  to  say  that  distributed  forcing  produces  a  strong  coupling  between  the 
induced  flow  and  the  forcing.  Perhaps  this  significant  result  has  been  appreciated  for  some 
time  by  experimentalists  who  have  induced  Gortler  vortices  by  placing  strips  of  tape  on 
the  wall  of  the  test  section  of  the  experiment. 

Now  let  us  consider  some  results  we  have  obtained  for  a  wall  with 

q{z)  =  ^e^.  (4.20) 

4 

In  order  to  monitor  the  strength  of  the  induced  vortex  activity  we  computed  u*{k,x,y)  the 
maximum  value  of  u  as  a  function  of  fc,  i  and  y  and  the  wall  shear  stress  as  a  function  of  k. 
By  inverting  the  transform  in  z  numerically  we  were  then  able  to  obtain  contours  of  ,  the 
inverse  of  u* ,  and  the  wall  shear  as  functions  of  x  and  z.  The  results  obtained  for  the  two 
functions  were  similar  so  we  shall  concentrate  on  results  for  uu  •  Figures  (4.6a,b,c,d)  show 
results  obtained  for  uu  for  the  obstacle  defined  by  (4.20).  The  results  shown  correspond 
to  G  =  4,  8,  12,  16.  The  size  of  G  is  a  measure  of  how  soon  we  expect  vortices  to  grow;  in 
the  absence  of  a  unique  neutral  curve  for  the  Gortler  problem  we  cannot  be  more  precise 
than  this.  In  the  horizontal  direction  x  varies  from  0.7  to  9  whilst  in  the  vertical  direction 
2  varies  between  -20,20.  We  could  have  demonstrated  the  same  effect  shown  in  Figure  (4.6) 
by  fixing  G  and  varying  the  position  where  the  forcing  was  switched  on.  Figure  (4.6)  shows 
clearly  that  immediately  after  the  obstacle  the  disturbance  velocity  field  decays  initially 
and  is  concentrated  in  a  wake  behind  the  obstacle.  Further  downstream  the  curvature  of 
the  wall  re-amplifies  this  disturbance  field  into  Gortler  vortices  which  spread  out  within  a 
wedge-shaped  region.  The  position  where  this  amplification  takes  place  is  a  function  of  G; 
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as  we  would  expect  amplification  occurs  closer  to  the  obstacle  as  G  is  increased.  We  note 
that  the  onset  of  vortices  in  such  wedge  shaped  regions  has  been  observed  expermentally  by, 
for  example,  Mangalam  et  al  (1987).  We  conclude  that  a  roughness  element  of  lengthscale 
0(L),  0(i2e"  a  X)  in  the  streamwise  and  spanwise  directions  respectively  will  produce  an 
induced  vortex  field  with  an  almost  negligible  drop  in  amplitude.  Moreover  the  onset  of 
vortex  activity  will  occur  in  a  wedge-shaped  region  behind  the  obstacle. 

Further  calculations  were  carried  out  for  asymmetric  obstacles  and  in  Figures 
(4.7a,b,c,d)  we  show  results  corresponding  to  Figure  (4.6)  but  with  q  now  given  by 


Q  = 


TT^ze  i’ 

Z2~ 


(4.21) 


The  conclusions  given  above  for  the  symmetric  case  clearly  also  apply  to  (4.4),  the  only 
difference  is  that  the  vortex  activity  is  now  symmetric  about  z  =  0. 


5.  The  most  unstable  Gortler  vortex 

We  recall  that  the  only  part  of  the  wavenumber- Gortler  number  plane  which  can  be 
described  in  a  self-consistent  manner  by  a  parallel  flow  theory  corresponds  to  large  values 
of  the  Gortler  number  and  wavenumber.  More  precisely  we  note  that,  in  the  limit  G  — »  oo 
with  A:  ~  G«,  neutral  vortices  described  by  the  asymptotic  structure  of  Hall  (1982a) 
exist.  Moreover  in  the  neighbourhood  of  the  neutral  state  spatially  amplified  modes  with 
growth  rates  0(A:^)  occur.  These  modes  develop  in  the  streamwise  direction  in  such  a 
way  as  to  maximise  their  local  growth  rate.  In  order  to  isolate  the  fastest  growing  mode 
at  high  Gortler  number  (note  that  this  concept  is  not  tenable  at  finite  Gortler  numbers) 
we  investigate  the  wavenumber  spectrum  close  to  the  right  hand  branch  of  the  neutral 
curve.  The  analysis  is  based  almost  entirely  on  that  given  by  Hall  (1982a),  so  we  give  only 
the  most  relevant  results  here;  in  fact  the  analysis  of  §3  essentially  outlines  the  neutral 
structure  of  such  modes  if  the  forcing  is  switched  off. 

Consider  then  the  solution  of  (2.12a),  (2.12g)  in  the  limit  G  —*  oo  with  k  =  0(G*). 
More  precisely  we  write 

k  =  AG‘  (5.1) 

where  A  =  0(1)  initially  but  later  it  will  be  informative  to  move  further  away  from  the 
right  hand  branch  by  letting  A  — »  0.  We  seek  a  solution  of  (2  12a),  (2.12g)  in  a  layer  of 
depth  fc"  >  ~  G~  •  centred  on  the  location 

y  =  V{x). 
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The  appropriate  expansion  of  the  disturbance  velocity  u  there  is 

u  =  e°  ^  ^  +  A:"  >  x)  +  •  •  •}  (5.2) 

together  with  a  simialr  expansion  for  k~\;  here  ^  =  A»  G>  {y  — and  we  have  anticipated 
a  growth  rate  of  size  0(Ga ).  In  the  layer  of  vortex  activity  we  write 

u  =  Uq{x)  +  \~ '^ni{x)ipG~  •  +  ■  ■  ■  (5.3) 


and,  following  Hall  (1982a),  we  find  that  the  zeroth  order  approximation  to  the  eigenrela- 
tion  is 

{Go/3  +  =  XQo(x)Gi(x).  (5.4) 

At  this  stage  y(x)  has  not  been  fixed  so  we  must  interpret  (5.4)  as  an  equation  to  determine 
f3  —  /3(A,x,  y(x)).  If  there  is  no  forcing  the  vortex  must  decay  exponentially  as  \ip\  — >  oo, 
this  provides  a  second  condition 


d 


^[(Q(x,y)/3  +  -  xHx,y)u^{x,y)]\  =  0. 

oy 


(5.5) 


The  vertical  structure  of  the  vortex  is  then  obtained  in  terms  of  parabohc  cylinder  functions 

as  was  found  in  §3.  We  can  then  eliminate  the  dependence  of  /?  on  y  by  combining  (5.4) 

and  (5.5)  to  give  /?  =  /?(A,x).  We  note  that  the  vortices  are  responding  in  a  quasi-parallel 

manner  so  we  could  in  fact  scale  x  out  of  the  problem  by  redefining  A  and  /?.  In  particular 

if  Q  corresponds  to  a  Blasius  boundary  layer  so  that  U  =  f'{r])  —  f'[  . — )  where  /  is  the 

V  2x 


Blasius  function  then  we  must  write 
the  simplified  problem  for  /?t(At)  is 

(/'/?t  +  At’)2  =  /7", 
2r^t(/'/?t  + At’)  =  (/'/")', 


in  which  case 


(5.6a, b) 


res 


where  the  Blasius  function  is  to  be  evaluated  at  some  77  =  yt  where  (5.6b)  holds.  In  Figu 
(5.1a,b)  we  show  j3^  and  ryt  as  functions  of  At.  The  cut-off  value  of  r^t  ~  1.5  corresponds  to 
the  neutral  state  found  by  Hall  (1982a).  However,  we  see  that  /?t  is  a  monotonic  decreasing 
function  of  At  and  in  fact  has  a  singularity  when  At  — ♦  0-j-.  An  analysis  of  (5.6a,b)  shows 
that  in  that  limit 

^t  xt  (5.6c) 
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so  that  ^mall  wavelength  vortices  become  progressively  more  unstable  as  — ; - >  0.  Thus 

we  reach  the  conclusion  that  the  quasi-parallel  approach  of  Hall  (1982a)  does  not  capture 
the  fastest  growing  mode  and  that  the  most  unstable  modes  are  probably  to  be  found  con¬ 
centrated  in  a  layer  approaching  the  wall.  Interestingly  it  is  found  that  the  fastest  growing 
temporal  mode  is  captured  by  the  asymptotic  structure  of  Hall  (1982a);  unfortunately  the 
temporal  modes  are  irrelevant  to  the  receptivity  problem  for  free-stream  or  wall  induced 
vortices. 

The  above  discussion  suggests  that  the  most  amplified  wave  for  G  >>  1  is  to  be  found 
where  k  <C  G  * .  Suppose  then  that  in  order  to  find  out  the  precise  location  of  this  mode  we 
investigate  the  region  G  1,A:  =  0(1).  In  fact  here  the  vortices  satisfy  the  inviscid  form 
of  the  vortex  equations  and  have  growth  rates  0(Ga).  We  are  unaware  of  any  previous 
work  on  the  spatial  modes  of  the  inviscid  Gortler  problem  so  we  give  slightly  more  detail 
now  than  was  given  for  the  fc  ~  G«  problem.  The  temporal  inviscid  Gortler  problem  is 
virtually  idential  to  the  inviscid  Taylor  problem  and  is  set  as  an  exercise  by  Drazin  and 
Reid  (1979),  and  discussed  by  Floryan  (1986).  Now  let  us  look  at  the  structure  of  the 
spatial  inviscid  Gortler  modes  and  see  if  they  include  the  fastest  growing  mode.  We  first 
expand  u  in  the  form 

u  =  6°"  ■^^‘'*{uo(>c.y)  +  G~^ui{x,y)  +  •••} 

together  with  a  similar  expansion  for  G~  >  no-  The  eigenvalue  problem  for  the  amplification 
rate  p  is  then  found  to  be 

-P‘^u[{d^  -  k'^)vQ  -  ^vol  = 

no  =  0,  y  =  0.  o~. 

Notice  that  since  (5.8)  corresponds  to  the  inviscid  limit  we  cannot  satisfy  the  viscous 
boundary  condition  no^  =  0.  For  a  given  value  of  x  we  can  solve  (5.8)  for  the  unstable 
eigenvalues,  corresponding  to  /?  >  0.  In  particular,  we  note  that  (5  8)  has  the  exact  solution 

VQ  -  Uexp{-ky),  p"^  =  kx/2,  (5.9a) 


valid  for  all  k. 

The  above  exact  solution  is  valid  only  if  Q  does  not  vanish  in  (0,oo)  and  corresponds 
to  an  unstable  mode  if  x  is  positive  so  that  the  wall  is  concave.  Our  primary  concern  in 
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this  paper  is  with  the  receptivity  problem  for  a  Blasius  boundary  layer  but  we  shall  now 
make  a  few  remarks  about  the  relevance  of  the  above  solution  to  other  flows. 

Firstly  we  note  that  a  second  exact  solution  of  the  diiferential  equation  for  is 

uo  =  nexp{ky),  =  —kx/2,  (5.96) 

but  this  solution  does  not  satisfy  the  required  boundary  condition  at  inflnity  unless  a 
vanishes  there.  More  precisely  we  require  that  Q  goes  to  zero  more  quickly  than  does  e*'' 

.  Thus  the  second  exact  solution  is  relevant  to  for  example  wall- jet  flows  and  corresponds 
to  an  unstable  disturbance  if  x  is  negative  so  that  the  wall  has  convex  curvature.  Hence  a 
wall-jet  flow  is  unstable  in  the  presence  of  either  convex  or  concave  curvature.  However  the 
magnitude  of  the  growth  rate  depends  only  on  the  magnitude  of  the  curvature,  whilst  the 
form  of  the  eigenfunction  depends  ony  on  the  sign  of  the  curvature.  We  further  note  that 
the  second  exact  solution  is  also  valid  only  for  flows  which  do  not  have  a  change  of  sign 
in  u  within  the  flowfield.  As  one  would  anticipate  the  unstable  modes  for  a  flow  having 
a  change  of  sign  of  U  are  more  complicated;  in  fact  they  have  a  critical  layer  structure 
where  u  vanishes  and  viscous  effects  must  be  retained  to  remove  the  singularity  in  the 
disturbance  velocity  field.  In  that  case  the  exact  solutions  describe  the  disturbance  only  in 
the  region  above  the  critical  layer  furthest  from  the  wall.  We  note  here  that  the  eigenvalue 
of  the  unstable  mode  associated  with  such  flows  is  still  as  given  by  (5.9a,b)  dependent  on 
the  sign  of  the  curvature.  A  result  which  is  even  more  surprising  is  that  (5.9a,b)  apply  to 
compressible  boundary  layers  for  which  u  does  not  change  sign.  Hence  for  a  compressible 
boundary  layer  the  temperature  profile  has  no  direct  effect  on  the  growth  rate  of  what  is 
in  fact  the  most  unstable  inviscid  mode.  In  fact  the  exact  solution  is  only  the  first  mode  in 
an  infinite  hierarchy  of  modes  which  yield  unstable  eigenvalues  if  Uy  is  anywhere  positive. 
However,  we  were  unable  to  determine  exact  solutions  of  (5.8)  for  these  higher  modes. 
In  Figure  (5.2)  we  present  the  first  two  unstable  eigenvalues  of  (5.8);  we  have  restricted 
attention  to  the  Blasius  flow  and  have,  for  convenience,  chosen  x  =  ^.  From  (5.8)  we 
expect  that  for  small  k  the  higher  modes  have  unstable  eigenvalues  /?  ~  fc.  This  prediction 
is  born  out  by  the  results  shown  in  Figure  (5.2)  for  the  second  mode. 

We  see  from  Figure  (5.2)  that  there  is  not  a  fastest  growing  mode  predicted  by  our 
calculations  since  for  the  two  modes  calculated  /?  tends  to  infinity  when  A:  — >  oo.  In  fact 
it  is  crucial  to  find  this  large  k  structure  since  it  provides  the  vital  clue  in  the  search 
for  the  most  unstable  mode.  The  calculations  which  we  performed  in  order  to  construct 
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Figure  (5.2)  showed  that  the  eigenfunction  of  any  mode  becomes  concentrated  at  the  wall 

d 

as  A:  — +  oo.  This  suggests  that  a  new  structure  will  emerge  from  (5.8)  when  —  ~  A:  in 

oy 

which  case  there  is  a  wall  layer  of  depth  A:“^  near  y  =  0.  Since  U  y  for  y  small  it  is  clear 
that  such  a  structure  can  exist  only  if  ~  A:> .  We  therefore  write 

ij)  =  ky,  /?  =  /?^A:> 


and  the  reduced  problem  for  /jt  is  then 


no 


}  = 


UQ  =  0,  x/)  =  0,  oo. 

This  is  a  form  of  Whittakers  equation 

W"  +  [-\  +  —  +  W  =  0 


(5.10) 


The  solutions  of  this  equation  using  the  notation  of  Abramowitz  and  Stegun  (1965)  are 
^K,n  (•2^))  (2)  and  can  be  expressed  in  terms  of  Kummer  functions.  The  above  equation 


reduces  to  (5.10)  if  /i  =  ±5  and 


2/?t2 


X 
K  • 


We  can  infer  from  the  small  and  large  z  forms  of  M,  W  given  by  Abramowitz  and 
Stegun  that  if  we  are  restricted  to  y,  =  then  the  only  values  of  K  which  enable  W  to 
vanish  at  2  =  0,  00  are  K  =  1, 2,  3  . . ..  Thus  the  unstable  values  of  /jt  are  given  by 

V!’  V5-- 

Note  here  that  the  first  eigenvalue  is  the  exact  solution  valid  for  all  k.  Thus  in  the  high 
wavenumber  limit  the  inviscid  modes  become  progressively  more  unstable  with  growth  rate 
~  G  3  A:  j .  On  the  other  hand  we  see  that  the  viscous  right  hand  branch  modes  become  more 
unstable  as  k  decreases,  actually  from  (5.1),  (5.2)  and  (5.6c)  we  see  that  in  that  regime  the 

I  k  2 

growth  rate  is  ~  G^i — j-)”  .  This  suggests  that  the  viscous  and  inviscid  modes  in  these 

G^ 

limiting  cases  will  merge  when  k  ~  G»  since  they  will  both  then  predict  growth  orders 
~  G* .  Moreover,  both  limiting  forms  then  suggest  that  the  vortices  will  then  be  confined 
to  a  layer  of  thickness  G»  at  the  wall.  Thus  in  this  ‘overlap  regime’  we  write 

k  =  AG^ 

=  ky, 
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and  expand  u  in  the  form 


3 

together  with  a  similar  expansion  for  G«vo-  Near  the  wall  the  basic  flow  Q  ~  p(x)y  and 
the  function  pair  {uq^vq)  is  found  to  satisfy 


r  d2  A  __ 

ld^2  x3 


ld^2  p 

/  d2  A!  A  XM^UO 

ld^2  p  nd^2  X3 

UQ  =  VQ  =  Vq  =  0,  '^  =  0,  OO. 


(5.11) 


The  limiting  inviscid  form  is  then  obtained  from  (5.11)  by  taking  the  further  limit  A  — »  0 
in  (5.11),  whilst  the  limiting  form  (5.6)  is  recovered  by  allowing  A  — >  oo  in  (5.11).  Thus 
without  further  calculation  we  know  that  the  eigenvalues  of  (5.11)  are  such  that 


/3~A>,  A-»0,  /3~A-2,  A-^oo. 


(5.12) 


so  that  at  some  intermediate  value  of  A  a  most  amplified  case  will  exist.  In  order  to  confirm 
this  we  solved  (5.11)  numerically  to  obtain  the  first  few  unstable  eigenvalues  as  functions 
of  A.  For  computational  purposes  it  is  convenient  to  eliminate  the  parameters  p,x  from 
(5.11)  by  writing 


A  =  (xm'^)‘A,  f3  =  {x^n)>P,  UQ  =  nuo,  uo=(x/^^)‘vo 


(5.13a, b,  c,  d) 


in  which  case  (5.11)  may  be  rewritten  as 


r  d2  _  VO 

UQ  =  5o  =  Vq  =  0,  ’5’  =  0,  OO. 


'Fuo 

T^’ 


(5.14) 


Thus  (5.14)  now  constitutes  an  eigenvalue  problem  0  =  /9(A)  and  the  first  two  unstable 
eigenvalues  are  shown  in  Figure  (5.3);  we  see  that  (5.11)  is  confirmed  and  that  each  unstable 
mode  has  its  growth  rate  maximised  at  a  finite  value  of  A.  Our  calculations  showed  that 
the  most  unstable  mode  corresponds  to  A  =  Ac  =  A76,P  =  .312.  Actually,  in  addition  to 
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the  modes  shown  there  is  an  infinite  sequence  of  less  unstable  modes  with  a  similar  P  —  X 
structure  moved  progressively  closer  to  the  A  axis. 

Thus  the  above  calculation  has  shown  that  the  most  unstable  linear  Gortler  vortex  at 
high  Gortler  numbers  is  viscous  but  with  wavenumber  ~  G»  rather  than  ~  G*  as  would  be 
appropriate  to  the  unstable  modes  close  to  the  neutral  curve.  A  further  significant  result  is 
that  the  most  unstable  vortices  are  close  to  the  wall;  this  leads  us  to  expect  that,  when  the 
appropriate  receptivity  problem  for  these  modes  is  discussed,  0(1)  ‘coupling  coefficients’  are 
possible.  This  problem  is  addressed  in  the  next  section;  to  close  this  section  we  make  a  few 
remarks  about  the  nature  of  the  eigenfunctions  associated  with  the  most  dangerous  mode 
we  have  isolated.  The  uq,  vq  velocity  fields  of  the  most  dangerous  mode  are  shown  in  Figure 
(5.4).  We  see  that  they  have  the  characteristic  shape  of  vortices  measured  experimentally. 
In  fact  Figure  (5.5)  shows  a  comparison  of  this  most  unstable  eigenfunction  with  that 
measured  by  Swearingen  and  Blackwelder  (1987).  Also  shown  are  the  eigenfunctions  from 
Smith’s  (1955)  parallel  flow  theory  and  a  nonparallel  solution  obtained  by  Hall  (1983).  We 
see  that  the  predicted  velocity  field  is  close  to  that  measured  experimentally. 

Now  let  us  consider  the  evolution  of  a  fixed  wavelength  vortex  as  it  develops  in  the  x 
direction.  We  note  that  it  is  observed  experimentally  that  Gortler  vortices  develop  such 
that  their  wavelength  is  conserved  so  A  =  constant  describes  the  experimental  situation. 
For  definiteness  we  suppose  that  the  basic  flow  is  a  Blasius  boundary  layer  so  that 

.  _i  .4696  _i 


and  further  that  the  wall  shape  is  such  that 

X  ~  ,  with  N  fixed. 

We  can  determine  /?(A)  for  such  a  situation  by  noting  that 

P  =  jj)  '  P 


Experimentally  it  is  observed  that  when  a  vortex  develops  downstream  it  conserves  its 
wavelength.  Thus  for  a  given  fixed  value  of  N  we  compute  A  =  A(x  "~g^ )  ‘  a^t  the 

location  x  and  read  off  the  corresponding  value  of  p  from  Figure  (5.3).  The  growth  rate 
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/3  at  that  location  is  then  computed  from  (3  =  ‘/3.  It  follows  that  the  maximum 

possible  growth  rate  at  any  x  will  increase  with  x  only  if  iV  >  ^.  However  we  note  that, 
in  the  special  case  =  1,A  is  a  constant  so  that  the  most  unstable  mode  occurs  at  the 
same  value  of  A  for  all  x  and  has  magnitude  increasing  like  .  This  particular  choice  of 
N  =  1  is  significant  for  the  asymptotic  structure  we  have  developed  because  in  this  case 
the  local  Gortler  number  varies  like  the  fifth  power  of  the  local  wavenumber  so  that  only 
the  magnitude  of  the  maximum  growth  rate  can  depend  on  x.  We  note  in  passing  that 
=  2  plays  a  similar  role  for  neutral  vortices  corresponding  to  the  right  hand  branch  of 
the  neutral  curve. 

In  fact  if  we  do  not  take  A^  =  1  we  see  that  A  either  decreases  or  increases  with  x 
depending  on  whether  N  >  1  oi  N  <  1.  This  means  that  for  a  given  value  of  A  the  growth 
rate  P{x)  will  have  a  maximum  value  at  some  value  of  x  if  iV  ^  1.  Thus  if  A/^  ^  1  there  is 
a  downstream  location  where  a  vortex  of  given  wavelength  is  most  amplified.  In  Figures 

o 

(5.6a,b,c)  we  show  /?  as  a  function  of  A  and  x  for  the  cases  N  =  0,-^,  1.  It  should  be  noted 
that  the  development  of  a  given  vortex  can  be  found  from  these  figures  by  taking  a  slice 
through  that  figure  parallel  to  the  (3  —  x  plane. 


6.  The  receptivity  problem  for  the  most  unstable  Gortler  mode 

We  have  seen  above  that  the  crucial  property  of  the  most  unstable  vortex  is  that  it 
is  confined  to  a  region  of  depth  G~  •  near  the  wall.  Suppose  then  that  the  wall  forcing 
function  F{x)  is  of  the  form 


F  =  0,  X  <  X, 

F  =  F{JG»[x  —  x]),  X  >  X. 


(6.1) 


Here  we  have  introduced  the  factor  J  —  (xV)*  for  convenience.  Thus  we  are  now  allowing 
the  wall  forcing  function  to  vary  on  a  length  scale  comparable  to  the  length  scale  over  which 
the  most  unstable  streamwise  vortex  develops.  If  we  then  define  x  by 


X  =  JG*  {x  —  x) 


and  write 

u  =  uo(x,  "F)  +  G“  •  ui(x,  ^)  +  •  •  • , 
{xfJ.G)~'nv  =  vo{x,^)  +  G~  •vi{x,<i!)  +  ■  ■  • , 


(6.2) 
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then  it  follows  directly  from  (5.14)  that  uq.vo  satisfy 


{ 

{ 
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1  a  I  1 


(6.3a, 6) 


Here  A,  are  as  defined  in  the  previous  section, 
conditions 

UQ  -  F(x),  VQ  =  0, 
— *  0.  ^  »  00. 


The  forced  problem  then  satisfies  the 
^'  =  0, 

(6.4a, b) 


In  fact  we  could  modify  the  boundary  condition  as  ^  »  00  in  order  to  allow  for  the 

possibility  that  vortices  are  induced  by  free-stream  disturbances.  However,  since  the  dis¬ 
turbance  in  the  free-stream  must  communicate  with  the  wall  layer  over  an  0(1)  distance  on 
a  G“  •  lengthscale  through  a  WKB  structure  similar  to  that  of  §3  the  coupling  coefficient 
for  such  a  disturbance  would  be  0(6"^  * )  and  so  this  type  of  excitation  will  be  weak  com¬ 
pared  to  that  discussed  here.  Alternatively  we  could  investigate  the  effect  of  wall  suction 
by  setting  uq  =  0  at  the  wall  and  applying  a  nonzero  boundary  condition  on  uo  there.  The 
forced  problem  (6.3)-(6.4)  is  most  easily  solved  by  again  taking  a  Laplace  transform  with 
respect  to  the  streamwise  co-ordinate.  If  we  denote  the  transforms  of  uo.uOt  and  F  by 
ao,Vo  and  F,  and  assume  that  there  is  no  incoming  vortex  field  then  we  find  that  uq,  uo 
satisfy  (5.14)  with  /3  replaced  by  the  transform  variable  s  and  the  boundary  conditions 
unchanged  except  for  replacing  the  condition  Go  =  0,  =  0  by  Go  =  F(s),  ^  =  0.  We  can 

write  the  solution  of  this  problem  in  the  form 


(QO,Go)  =  ^(5)(Ti*(s,4'),U-(s,«-)) 


(6.5) 


where  (u*,v')  satisfies  the  same  problem  as  (Go,  Go)  but  with  F  =  1.  It  is  obvious  at  this 
stage  that  u*  and  v'  will  have  singularities  in  the  complex  plane  where  3  =  0,,  the 

eigenvalue  of  (5.14)  for  j  =  1,2 _  Moreover  when  we  invert  the  transformed  velocity 

field  these  singularities  will  in  effect  fix  the  form  of  the  induced  vortex  field.  It  is  therefore 
essential  that  we  should  find  the  form  of  these  singularities.  We  can  obtain  the  form  of 
{u*  ,v')  near  the  singularity  by  expanding  {u*,v*)  in  the  form 


(^•,t^*) 


(^0>  ^0) 

3  -  0j 


+  uj)  +  ■  ■  •  . 


(6.6) 
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If  the  above  expansions  are  substituted  into  the  differential  system  for  (u*,u*)  and  s  is 
replaced  by 

5  =  +  (6.7) 

and  terms  of  order  {s  —  ^j)~^  are  equated  in  that  system  we  obtain 

(uq.vq)  =  Aj{uQj0j,y),voj0j,y))  (6.8) 


where  A,  is  a  constant  and  (u0;,U0j)  is  the  eigenfunction  of  (5.14)  appropriate  to  /3  = 

At  next  order  in  the  expansion  procedure  we  find  an  inhomogeneous  form  of  (5.14)  to 
determine  (ui,ui).  If  we  define  (by,/,)  to  be  the  function  pair  adjoint  to  (uoj,t^Oj)  then 
{bj ,  fj )  satisfies 


JL 

^3-  ^2’ 

bj  =  fj  =  /,-  =  0.  ^  =  0.  oo. 

and  the  solvability  condition  for  the  system  for  (ui,Ui)  yields 


A3  ■  A3 


(6.9) 


Ay  = 


J"  ^  ^  i  ^(  ~  t)*  0  J 


X> 


+ 


4’u.j 

A» 


(6.10) 


For  convenience  we  assume  at  this  stage  that  uoy  has  been  normalized  to  have  a  maximum 
value  of  unity  so  that  Ay  is  a  measure  of  the  transformed  disturbance  amphtude.  It 
remains  now  for  us  to  invert  the  transformed  velocity  field;  in  order  to  do  so  we  must 
be  more  precise  about  the  function  F.  If  we  are  interested  in  the  effect  of  an  isolated 
roughness  element  we  can  take  F  to  be  a  function  of  compact  support  and  then  F  cannot 
have  any  singularities  in  5,  >0.  If  we  wish  to  consider  the  effect  of  distributed  roughness 
this  restriction  must  be  relaxed  and  we  merely  s^sume  that  F  exists.  In  the  first  instance 
we  assume  that  F  has  compact  support  then,  since  the  only  singularities  of  (5.10)  in  5,  >  0 
correspond  to  the  simple  poles  discussed  above  the  contour  of  integration  for  the  inverse 
transform  is  chosen  parallel  to  the  imaginary  axis  to  the  right  of  s  =  /?i.  The  contour  is 
then  closed  in  the  left  hand  half  place  Realms)  <  and  the  only  contributions  to  this 
integral  come  from  the  poles  s  =  Pj.  Thus  we  obtain 


(uQ,  uo)  =  E  (woy ,  voy  '  Ay  F0j ) 

y=i 


(6.11) 
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so  that  for  a  given  value  of  A  the  effective  coefficient  is  AjF{Pj)  for  the  amplitude  of  the 
mode. 

However,  it  is  clear  that  for  large  i  (6.11)  will  be  dominated  by  the  fastest  growing 
mode  which  of  course  corresponds  to  j  =  1.  Thus  when  the  Fourier  transform  is  inverted 
to  give  the  effect  of  a  particular  Z  structure  for  the  wall  forcing  by  deforming  the  contour 
appropriately  we  will  see  that  far  downstream  a  vortex  is  set  up  with  growth  rate  and 
with  the  wavelength  A  corresponding  to  0i.  At  this  stage  we  could  repeat  the  discussion  of 
§4  and  determine  the  effect  of  a  three-dimensional  obstacle  on  the  onset  of  Gortler  vortices. 
In  fact  the  development  of  vortices  in  this  region  is  essentially  identical  to  that  of  §4.  Thus 
behind  a  three-dimensional  obstacle  there  will  be  a  wake  where  initially  the  effect  of  the 
obstacle  is  concentrated;  further  downstream  this  wake  will  expand  into  a  wedge-shaped 
region  of  vortex  activity  with  a  vortex  wavelength  corresponding  to  the  most  unstable 
mode  of  (5.10).  The  above  discussion  has  assumed  that  F  has  compact  support;  if  that  is 
not  the  case  then  some  further  discussion  is  required.  In  fact  if  F  does  not  have  a  pole  to 
the  right  of  s  =  /3i  then  the  large  x  downstream  behaviour  of  the  vortex  flow  will  be  as 
described  above,  since  the  associated  velocity  field  will  still  be  dominated  by  the  pole  at 
s  =  /3i.  However,  if  F  has  a  pole  to  the  right  ois  =  0i  at,  say,  s  =  s*  then  F  increases  like 
<2*'*  for  large  x  and  our  original  assumption  that  the  wall  perturbation  is  small  becomes 
invalid.  Thus  the  only  situations  described  by  our  theory  have  the  inverted  velocity  field 
dominated  by  the  fastest  growing  mode  of  (5.10).  Hence  the  question  of  whether  the 
roughness  element  is  isolated  or  distributed  is  not  resolved  for  the  fastest  growing  mode  in 
either  case.  In  either  case  F0i)  will  be  an  0(1)  number  so  that  each  type  of  forcing  leads 
to  an  0(1)  coupling  coefficient  between  the  surface  perturbation  and  the  induced  vortex 
field.  In  Figure  (6.1)  we  present  a  plot  of  Ai  versus  A;  we  can  see  from  this  plot  that  Ai 
attains  its  maximum  at  a  finite  value  of  A.  However,  there  is  no  reason  to  associate  this 
value  of  A  at  which  the  maximum  of  Ai  occurs  with  that  found  in  §5  for  the  most  unstable 
vortex. 

Finally  in  this  situation  we  make  a  few  remarks  about  more  similarities  between  the 
solutions  of  (6.3)  and  those  obtained  in  §4.  Firstly  we  note  that  (6.3)  has  a  small  x 
similarity  solution  essentially  identical  to  that  found  in  §4.  Again  the  similarity  variable 
is  C  =  -^  and  for  ^  =  0(x>)  it  is  easy  to  show  that  the  solution  expands  as 

X  a 

uo  =  rGo(C)  +  •  •  ■ ;  VO  = -:Vqo(C)  +  •  •  •  • 

X  x» 
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with  (tio,  given  from  a  rescaled  form  of  (4.3).  Above  this  layer,  where  =  0(1),  there  is 
a  region  with  uq  ~  ■:T^o(y),vo  ~  -T®o(y)  again;  in  fact  it  turns  out  that  So  ~  e"**'  there 

X  i  S  3 

so  that,  if  the  2  structure  of  the  wall  perturbation  is  concentrated  near  2  =  0  and  modelled 

by  a  delta  function,  above  the  layer  the  normal  perturbation  velocity  in  x,y,z  space 

is  proportional  to  x~»—k - o-  Note  that  in  the  wall  layer  the  spanwise  structure  of  the 

disturbance  remains  concentrated  in  a  delta  function  form.  However  the  similarity  form 
z 

~  —55 - s  shows  clearly  the  ‘wake’-structure  induced  by  a  localized  disturbance.  This  is 

precisely  the  situation  observed  in  the  numerical  investigations  of  §4. 

7.  Conclusions 

We  shall,  in  the  first  place,  draw  some  conclusions  about  the  stimulation  of  Gortler 
vortices  by  wall  roughness  with  spanwise  wavelength  comparable  to  the  boundary  layer 
thickness.  The  key  point  in  this  situation  is  that  Gortler  vortices  develop  in  a  non-parallel 
manner;  there  is  no  such  thing  as  a  unique  growth  rate  for  a  growing  vortex  or  a  unique 
neutral  curve  for  a  vortex  which  is  not  growing  or  decaying.  Necessarily,  this  means  that 
the  receptivity  problem  for  vortices  of  0(1)  wavelength  must  be  dominated  by  non-parallel 
effects.  In  §4  we  discussed  the  0(1)  wavelength  receptivity  problem  for  roughness  elements 
varying  in  the  streamwise  direction  on  a  lengthscale  comparable,  and  shorter  than,  the 
body  scale.  However  in  the  latter  case  our  analysis  fails  if  the  roughness  varies  in  the 
streamwise  direction  on  the  boundary  layer  lengthscale.  We  saw  in  §4  that,  for  a  forcing 
function  of  the  form  (4.1)  with  F  being  a  function  with  compact  support,  the  effect  of 
the  forcing  is  to  generate  a  similarity  solution  of  the  linear  Gortler  vortex  equations  in  the 
neighbourhood  of  where  the  forcing  occurs.  By  integrating  the  disturbance  equations  with 
the  similarity  solution  used  to  generate  the  initial  vortex  field  we  were  able  to  determine 
‘the’  neutral  curve  appropriate  to  an  isolated  roughness  element  inserted  into  the  boundary 
layer.  Furthermore  we  found  that  this  type  of  forcing  is  a  relatively  inefficient  generator 
of  Gortler  vortices  because  the  vortex  amplitude  at  the  position  where  it  begins  to  grow 
is  an  order  of  magnitude  smaller  than  the  amplitude  of  the  forcing. 

Where  the  forcing  varies  on  the  body  lengthscale,  the  coupling  coefficient  is  0(1).  In 
this  case  we  found  that  behind  a  roughness  element  localized  in  the  spanwise  direction  the 
disturbance  initially  decayed.  Further  downstream  the  disturbance  was  found  to  be  ‘re¬ 
amplified’  by  centrifugal  effects  and  a  wedge  of  vortex  activity  generated.  This  pattern  of 
vortex  activity  has  been  observed  experimentally  by,  for  example,  Mangalam  et  al  (1987). 
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In  §3,  §5  and  §6  we  concentrated  on  the  case  where  the  Gortler  number  is  large  and 
the  vortex  wavelength  is  small.  The  investigation  of  §4  showed  conclusively  that  small 
wavelength  vortices  associated  with  the  right  hand  branch  of  the  neutral  curve  suffer 
an  exponential  drop  in  amplitude  because  the  region  where  they  are  unstable  is  an  0(1) 
distance  from  the  region  where  the  forcing  occurs.  Almost  certainly  this  means  that  it  is 
unlikely  that  small  wavelength  Gortler  vortices  could  be  produced  by  wall  roughness  in  an 
experiment. 

However,  there  is  a  situation  when  small  wavelength  Gortler  vortices  will  be  preferen¬ 
tially  induced  by  wall  roughness.  We  refer  to  the  region  discussed  in  §5,  6,  where  the  vortex 
wavelength  associated  with  the  forcing  is  0(G» ).  It  is  interesting  to  note  that  the  fastest 
growing  Gortler  vortex  described  in  §6  had  not  been  identified  by  previous  authors  because 
it  had  been  assumed  that  such  a  mode  would  have  to  be  associated  with  disturbances  close 
to  the  right-hand  branch  of  the  neutral  curve.  In  fact,  the  fastest  growing  mode  occurs 
in  an  overlap  region  between  inviscid  modes  of  wavenumber  0(1)  and  right-hand  branch 
modes  with  wavenumber  0(^4).  The  eigenvalue  problem  (5.11),  which  determines  the 
growth  rate  of  the  fastest  "rowing  mode,  is,  in  fact,  clearly  related  to  the  Taylor  vortex 
problem  for  the  flow  0*=’  veen  concentric  cylinders  with  the  outer  cylinder  at  rest  and  the 
inner  cylinder  rotf'.ting.  In  that  case  corresponds  to  the  azimuthal  amplification  rate  of 
a  steady,  smsli  wavelength  Taylor  vortex  trapped  near  the  outer  cylinder.  However  (5.11) 
has  no  neutral  modes  so  that,  for  the  Taylor  problem,  where  the  disturbance  must  be 
single- valued  in  the  azimuthal  direction,  this  type  of  disturbance  is  not  physically  relevant. 
However,  in  other  situations  where  no  such  constraint  is  appropriate  (5.11)  determines  the 
fastest  growing  spatial  mode  of  instability.  Thus  the  eigenvalue  problem  (5.11)  is  directly 
relevant  to  the  Taylor-Dean  instability  in  pressure  gradient  driven  flows  in  curved  geome¬ 
tries,  whilst  a  generalization  of  (5.11)  involving  a  further  second  order  equation  is  relevant 
to  the  compressible  Gortler  problem  with  and  without  curvature.  We  see  from  (5.11)  that 
the  only  dependence  of  the  amplification  rate  on  the  mean  state  is  through  the  wall  shear 
fi,  and  in  fact  /?  ~  ^».  Thus  the  stabilizing  effects  of,  for  example,  wall  suction  can  be 
inferred  directly  from  their  effect  on  the  wall  shear. 

Finally,  we  turn  to  the  results  we  obtained  in  §6  for  the  receptivity  problem  associated 
with  the  fastest  growing  mode.  Our  main  result  in  that  section  was  that  Ai,  the  coupling 
co-efficient  of  the  fastest  mode,  is  an  0(1)  quantity.  This  means  that  wall  roughness  at  high 
Gortler  numbers  will  stimulate  vortices  of  wavenumber  G»  close  to  the  position  where  the 
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roughness  begins.  If  the  roughness  is  localized  in  the  spanwise  direction  the  vortex  activity 
again  begins  in  a  wedge  behind  the  obstacle.  It  remains  an  open  question  whether  finite 
amplitude  effects  will  significantly  influence  the  structures  we  have  discussed  in  this  paper. 
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Figure  (3.1)  The  different  regions  in  the  high  wavenumber  limit  receptivity  problem. 


Figure  (4.4a)  The  downstreaim  development  of  the  velocity  field  induced  by  the  similaunty 
solution  for  the  caise  k  =  .45,  G  =  S.  The  disturbance  corresponds  to  a  delta  function  forc¬ 


ing  fvmction  at  x  =  .5  ajid  is  shown  at  i  =  .521,  .541,  .561,  .581,  .601, 1.401, 2.201, 3.001, 3.801. 
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Figure  (4.4b)  The  downstream  development  of  the  velocity  field  induced  by  the  similarity 
solution  for  the  case  k  =  .45,  G  =  8.  The  disturbamce  corresponds  to  a  delta  function  forc¬ 


ing  fimction  at  I  =  .5  and  is  shown  at  i  =  .521,  .541,  .561,  .581,  .601, 1.401, 2.201, 3.001, 3.801 
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Figure  (4.6)  Contour  plots  of  in  the  x-z  plane  for  the  obstacle  defined  by  (4.20).  The 
results  correspond  to  a)  G  =  4,  b)  G  =  8,  c)  G  =  12  and  d)  G  =  16. 
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Tlic  clepeiideiice  of  the  the  first  two  unstable  eigenvalues  (i  of  (5.8)  on 
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.5)  A  coinj)arisoii  between  the  expeiiuientally  measured  eigenfunction  of  S 


le  'cuuplitiule’  Aj  as  a  fiiiu  tiou  of  A. 
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16  Abstract 


The  receptivity  probiem  for  Gortler  vortice*  induced  by  wall  roughness  is  investigated. 
The  roughness  is  modelled  by  small  amplitude  perturbations  to  the  curved  wall  over  which 
the  flow  takes  place.  The  amplitude  of  these  perturbations  is  taken  to  be  sufficient  y  smaii  I  ■ 
the  induced  Gortler  vortices  to  be  described  by  linear  theory.  The  roughness  is  issumed  to 
vary  in  the  spanwise  direction  on  the  boundary  layer  lengthscale,  whilst  in  the  flow  direction 
the  corresponding  variation  is  on  the  lengthsc^e  over  which  the  wall  curvature  varies.  In  fact 
the  latter  condition  can  be  relaxed  to  allow  for  a  faster  streamwise  roughness  variation  so  long 
as  the  variation  does  not  become  as  fast  as  that  in  the  spanwise  direction.  The  function  which 
describes  the  roughness  is  assumed  to  be  such  that  its  spanwise  and  streamwise  dependences 
can  be  separated;  this  enables  us  to  make  progress  by  taking  Fourier  or  Laplace  transforms 
where  appropriate.  The  cases  of  isolated  and  distributed  roughness  elements  are  investigated 
and  the  coupling  coefficient  which  relates  the  amplitude  of  the  forcing  and  the  induced 
vortex  amplitude  is  found  asymptotically  in  the  small  wavelength  limit.  It  is  shown  that 
this  coefficient  is  exponentially  small  in  the  latter  limit  so  that  it  is  unlikely  that  this  mode 
can  be  stimulated  directly  by  wall  roughness.  The  situation  at  0(1)  wavelengths  is  quite 
different  and  this  is  investigated  numerically  for  different  forcing  functions.  It  is  found  that 
an  isolated  roughness  element  induces  a  vortex  field  which  grows  within  a  wedge  at  a  finite 
distance  downstream  of  the  element.  However,  immediately  downstream  of  the  obstacle  the 
disturbed  flow  produced  by  the  element  decays  in  amplitude.  The  receptivity  problem  at 
larger  Gortler  numbers  appropriate  to  relatively  large  wall  curvature  is  discussed  in  detail. 
It  is  found  that  the  fastest  growing  linear  mode  of  the  Gortler  instability  equations  has 
wavenumber  proportional  to  the  one  fifth  power  of  the  Gortler  number.  The  mode  can  be 
related  to  ooth  inviscid  disturbances  and  the  disturbances  appropriate  to  tl  e  right  hand 
branch  of  neutral  curve  for  Gortler  vortices.  The  coupling  co-efficient  betv  een  this,  the 
fastest  growing  vortex,  and  the  forcing  function  is  found  in  closed  form. 
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